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ABSTRACT 

We present a self-consistent comparison of the rise times for low- and high-redshift 
Type la supernovae. Following previous studies, the early light curve is modeled using a 
law, which is then mated with a modified Leibundgut template light curve. The best- 
fit t"^ law is determined for ensemble samples of low- and high-redshift supernovae by 
fitting simultaneously for all light-curve parameters for all supernovae in each sample. 
Our method fully accounts for the non-negligible covariance amongst the light-curve 
fitting parameters, which previous analyses have neglected. Contrary to Riess et al. 
(1999a), we find fair to good agreement between the rise times of the low- and high- 
redshift Type la supernovae. The uncertainty in the rise time of the high-redshift Type 
la supernovae is presently quite large (roughly ±1.2 days statistical), making any search 
for evidence of evolution based on a comparison of rise times premature. Furthermore, 
systematic effects on rise-time determinations from the high-redshift observations, due 
to the form of the late-time light curve and the manner in which the light curves of 
these supernovae were sampled, can bias the high-redshift rise-time determinations by 
up to I'^'g days under extreme situations. The peak brightnesses — used for cosmology 
— do not suffer any significant bias, nor any significant increase in uncertainty. 

Subject headings: supernovae: general — cosmology: observations 

1. Introduction 

Two independent research groups have presented compelling evidence for an accelerating uni- 
verse from the observation of high-redshift Type la supernovae (SNe la) (Perlmutter et al. 1999; 
Riess et al. 1998). These findings have such important ramifications for cosmology that every effort 
must be made to thoroughly test the calibrated standard candles on which they are based. Indeed, 
these groups, and others, are pursuing additional observations at both high- and low-redshift to 
confirm these results. There are programs in place aimed at reducing the statistical errors, testing 
systematic errors, limiting the amount of absorption due to grey dust (Aguirre 1999), and searching 
for signs of evolution as a function of redshift in SNe la. 

Recently Riess et al. (1999a) attempted to examine the question of whether the rise times of 
SN la evolve. They used new low-redshift SNe la light-curve photometry from Riess et al. (1999b) 
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to compare the mean rise time of these SNe la to a prehminary rise time for high-redshift SNe la 
given in a conference abstract by Groom (1998) and based on a composite Ught curve derived from 
Supernova Cosmology Project (SCP) observations. Riess et al. noted a 5.8-cj difference between 
the rise times from the low-redshift data and from the Groom (1998) prehminary analysis of high- 
redshift data, with the high-redshift supernovae having shorter rise times by 2.4 days. Based on 
this result, they suggested the possibility that SNe la undergo sufficient evolution to account for 
what has been interpreted as evidence for an accelerating universe. 

In what follows, we address major shortcomings of these earlier analyses which fundamentally 
alter the conclusion of Riess et al. (1999a). Specifically, the analysis method used in Groom (1998) 
to produce a high-redshift rise-time estimate is very different than that used to produce the low- 
redshift rise-time estimate of Riess et al. (1999a). Furthermore, both analyses neglected correlated 
uncertainties in the light-curve fit parameters, and amongst the light-curve data points, so neither 
of these analyses is complete. We also examine, in §3, the role of light-curve sampling differences 
between the low-redshift and high-redshift SN la observations and how they can conspire with 
systematic deviations from the fitted reference template — seen for normal SNe la — to shift the 
inferred rise time. In §4 we briefly discuss the (small) impact on the cosmological application of 
SNe la resulting from light-curve variations. We conclude in §5 with a summary of our results and 
a discussion intended to help guide future work on the question of whether SNe la evolve. 



2. Statistical Analysis of SNe la Rise Times 

2.1. Description of the Problem 

Figure 1 illustrates the full SN la template, ^^{t)^ normally used by the SCP, which is a 
modified version of the Leibundgut template (Leibundgut 1988; Perlmutter et al. 1997). The light- 
curve fitting parameters are the peak fiux, fmax, time of maximum, tmax, and light curve stretch, 
s. (Note that all time dependent quantities refer to the rest frame of the supernova.) Goldhaber 
(1998) has demonstrated the remarkable fact that the stretch method applies to the rising portion 
of SN la light curves as well as it applies to the declining portion (up to +25 days after maximum) 
to better than 2% of the peak flux. This has been confirmed for nearby SNe la by Riess et al. 
(1999b). One can represent the flux light curve, f{t), as follows: 

f{t) = fmax1p{{t - tmax)/s) 

This approach works well in the U—, B— and F— bands over the range —20 days < t — tmax < 
+25 days (see both Perlmutter et al. (1999) and Perlmutter et al. (1997) for a full explanation of 
the use of this approach). 

A meaningful comparison of rise times for low- and high-redshift supernovae requires that 
both datasets be fit with the same template, and that the fits be performed in a manner which 
fully accounts for the covariance between the light-curve fitting parameters and the calculated rise 
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Fig. 1.— The modified Leibundgut template used by the SCP (Perlmutter et al. 1999, 1997), 
with the standard hght-curve fitting parameters, fmax, tmax-, and s labeled. Also shown arc two 
examples from the grid of laws mated to the modified Leibundgut template used to perform the 
statistical fits discussed in the text. The law is parameterized using an explosion day, t^xp, and 
a date, tjom, at which the law mates with the modified Leibundgut template. 
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time. For the high-redshift data, accounting for covariance in the hght-curve fitting parameters 
is especially important since the uncertainties on individual data points are relatively large. Such 
uncertainties allow the fitted date of maximum light, tmaxi the peak brightness, fmax, and the 
light-curve width, s, to be changed in compensating ways to yield similarly good fits. Thus, these 
parameters are correlated, and since determination of the rise time or explosion date, texpi involves 
both s and tmax, it is incorrect to fit for these parameters while holding s and tmax fixed. 

Take for example the case where the fitted value of tmax, t'max^ too early by 1 day. The 
fitted value of s, s', will suffer a compensating increase by roughly 1/15 in an effort to fit the data 
on the fast-declining, well-sampled portion of the light curve at -|-10 <t — tmax < +20 days. The 
effective stretch-corrected epoch, tg = {t — tmax)/s-i of a point nominally at f — tmax — 

-20 days 

and for s = 1 would be incorrect by: 

Ai. = {t- tmax)/s -it- O/s' 

-20 -19 
= LOO - 107 = ^'^y^- 
Likewise, if t^^^. were 1 day after tmax, s' would be smaller than the true s, changing Atg by roughly 
+2.5 days. This is the principal mechanism by which uncertainties in the light-curve fit parameters 
propagate into increased uncertainty in SN la rise times. (Our Monte Carlo simulations in §3 bear 
this out.) If the uncertainties in tmax and s had simply been propagated as if they were independent, 
the assigned uncertainty would be 1.7 days, and the correlated nature of the uncertainties would 
be lost. It is true that a point at t — tmax = —20 days may also play some role in constraining s. 
However, for the datasets considered here the observations on the rising portion of the light curves 
are generally much less certain than those on the declining portion. 

The analyses presented in Groom (1998) and Goldhaber (1998) were designed to test the 
efficacy of the stretch technique when applied to the rising portion of the light curves of high- 
redshift supernovae, and to attempt to improve that portion of the SCP light curve template. The 
high-redshift data from the SCP were aligned to stretch-corrected epochs, ts, using tmax and s 
for each supernova determined from individual light-curve fits without exclusion of data from any 
light-curve epoch. Then a t^ rise-time model was fit to the ensemble pre-max data, with the final 
result quoted for t^ fits covering rest-frame epochs —21 to —10 days with respect to tmax- None 
of the uncertainty due to the light-curve fitting parameters was propagated into the final quoted 
rise-time uncertainty (Goldhaber 1999, private communication). The resulting fit was then used 
to develop a revised template, and the individual SNe la light curves were then re-fit to this revised 
template. 

Riess et al. (1999a) analyzed the low-redshift data very differently: they aligned their low- 
redshift data using tmax and s for each supernova as in the preliminary high-redshift analysis, but 
they used only data from —10 to -|-35 days to fit the light curves. After aligning the light curves, 
a t^ rise-time model was fit to the ensemble pre-max data, with the final result quoted for t^ fits 
covering rest-frame epochs —23 to —10 days. Following Riess et al. (1999b), the uncertainty in 
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imax and s was accounted for in Riess et al. (1999a) by increasing the uncertainties on the stretch- 
corrected light-curve photometry points. The modest contribution due to correlated uncertainties 
was not included. 

Both of these studies fixed tmax and s for the individual SNe la before fitting the model 
from which explosion dates were inferred. They propagated the uncertainties in the light-curve 
fits parameters in an incomplete and approximate way. Since this approach does not allow each 
individual supernova's light-curve fit parameters, fmax, tmax, and s, to adjust to give the best 
fit as different rise times are tested, the uncertainties quoted in these studies are likely to be 
underestimates. In addition, since the two studies fit to different time intervals of data, a comparison 
of the central values may not be self-consistent. 

2.2. Fitting Method 

The most assumption- free means of accounting for how the uncertainties in the fits to individual 
SNe la light curves affect the value and uncertainty of the rise time is to explicitly test various rise 
times to see how well the SNe la are able to adjust to give fits of similar quality. This is more 
accurate than, and avoids difficulties associated with, attempting to propagate uncertainties based 
on the covariance matrix determined at the best-fit value when dealing with complex parameter 
probability spaces, such as those which occur for some SNe la light curves dealt with here. This 
approach requires that a family of templates with different rise times be defined and fit to the entire 
photometric dataset for each SN la. 

Unfortunately, at present, very little light-curve data are available for determining a suitable 
early-epoch template for a SN la. Therefore we have constructed a grid of templates consisting of 
models starting with zero flux at an explosion epoch, texp, and joined to the modifled Leibundgut 
template at epoch, tjoin- A model can be justified under the conditions of uniform expansion 
and constant effective temperature from simple physics (see also Arnett (1982)). Two examples 
from this family of t^-model, texp, tjoin templates are shown in Figure 1, with the epochs texp and 
tjoin labeled. These can be compared to the modified Leibundgut template, which is known to be 
a reasonable approximation to the light curves of many SNe la (with the timescale stretched or 
contracted) . 

Note that the use of t^xp, tjoin to describe the early-epoch light curve is simply a reparame- 
terization of the ce,t^^p models (i.e., f{t) = a{t — texpY) used in previous studies, with the added 
constraint of continuity where the a, t^^.^ model ends and the modified Leibundgut template begins. 
Riess et al. (1999a) did not impose a continuity constraint since the fitting to the early stretch- 
corrected light curve with an a,tl^p model was performed after (some portion) of the original light 
curve was fit with another template. Groom (1998) and Goldhaber (1998) have an implicit continu- 
ity constraint in that they mated their best-fit t^^jp model to the remainder of their light curve when 
constructing each new template. In this paper the fit for the rise time and the overall light-curve 
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parameters is performed simultaneously. 

An added benefit of our parameterization is that texp, tjoin are more nearly orthogonal than 
a,tl^p. This is because the already-established modified Leibundgut template provides a strong 
constraint on the amplitude of a a, tf,^.p model at the point it crosses the modified Leibundgut 
template, a simply adjusts itself to satisfy this constraint as t^xp is changed. (This leads to the 
narrow, but strongly tilted, confidence regions in Figure 1 of Riess ct al. (1999a)). 

The fitting method we use integrates the probability [P oc exp(— x^/2); see Eq. 28.22 in Ceolin 
ct al. (1998) ] over the parameters fmax, tmax, and s separately for each supernova, at each value 
of texp, tjoin- The fits are performed in flux (rather than magnitudes); this allows the use of non- 
detections, these being the principal source of early-epoch data for the Perlmutter et al. (1999) 
high-redshift supernovae. Two alternative methods are used to perform the integrations over fmax , 
tmaxi and s. In the first, the integral over fmax is performed analytically and the subsequent 
integration over tmax^ and s uses the adaptive integration algorithm of Berntsen et al. (1991). The 
second method uses a grid of Afmax = 0.01, Atmax = 0.1 days and As = 0.01 centered on the 
averages of the best fit values over texp for the high-redshift SNe la. To account for the tightly 
constrained parameters of the low-redshift SNe la a hybrid technique is used in which the integral 
over fmax is done analytically and a grid of Atmax = 0.01 days and As = 0.0005 is used to integrate 
over tmax and s. In each, the limits are chosen such that the probabilities are negligible at the 
boundaries. We find excellent agreement between each of these methods. 

The end product is a map of probability over texp, tjoin for each supernova. These probability 
maps are then multiplied, then renormalized, for an ensemble of supernovae, e.g., the high-redshift 
supernovae from Perlmutter et al. (1999) or the low-redshift supernovae of Riess et al. (1999a), to 
determine the joint probability distribution function over texp, tjoin, P{texp, tjoin)', or after normal- 
izing over texp for each tjoin, the conditional probability distribution function for texp given tjoin, 

exp\tjoin)- 



2.3. Supernova Light-Curve Samples 

What we will hereafter refer to as the "low-redshift SNe la" sample consists of SN 1990N, 
SN 1994D, SN 1996bo, SN 1996bv, SN 1996by, SN 1997bq, SN 1998aq, SN 1998bu, and SN 1998ef, 
for which early-epoch light-curve photometry transformed to i?-band from unfiltered CCD images 
has been reported by Riess et al. (1999b). The early-epoch photometry was supplemented with 
data from Lira et al. (1998); Patat et al. (1996); Meikle et al. (1996); Riess et al. (1999c); Suntzeff 
et al. (1999); Jha et al. (1999); Riess et al. (1999b) to produce full S-band light curves extending 
over peak and beyond. Riess et al. (1999b) reports four early-epoch light-curve points (one an 
upper limit) for SN 1998dh, however we were unable to include this supernova since the subsequent 
light-curve photometry was unavailable. 

What we will hereafter refer to as the "high-redshift SNe la" sample consists of the 30 SNe la 
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from Perlmutter et al. (1999) having redshift 0.35 < z < 0.65, with the exception of SN 1997aj-'^. 
As defined, this sample satisfies the requirements that at least 60% of the light in the i?-band 
comes from the rest-frame B-hand and that at least 60% of the rest-frame i?-band light is included 
in the R-hand. Redshift limits satisfying these conditions were determined using the -B-band and 
i2-band filter responses given in Bessel (1990), along with spectra of normal SNe la as a function of 
light-curve epoch constructed by Nugent et al. (2000). These restrictions allow comparison with the 
low-redshift S-band photometry of Ricss et al. (1999a) while minimizing the potential uncertainties 
inherent in making large cross-filter K-corrections. Even so, K-correction uncertainties will be 
present for those supernovae near these redshift limits, as well as at very early times where few 
spectra are available from which K-corrections can be calculated. Note that most of the other eleven 
supernovae from Perlmutter et al. (1999) have complete light curves only in rest-frame F-band or 
[/-band, and therefore are unsuitable for determination of the S-band light-curve parameters. Also 
note that not all 30 high-redshift SNe la have equal value in determining the rise time. Only those 
that were fortuitously caught on the rise in the reference images of the search run can constrain 
this region of the light curve. 



2.4. Results of the Statistical Analysis 

Templates were generated for —29.9 < t^xp < —10.1 days, in steps of 0.2 days, and for —20 < 
tjoin < —4 days, in 1 day steps. Fitted templates were required to have tgxp earlier than tjoin- 
Figure 2 presents the results of these fits; shown are the 1-, 2-, and 3-a confidence regions for the 
conditional probability, P{texp\tjoin), for the high-redshift SNe la sample. Also shown are points 
which mark the most probable value of texp at each tjoin for the low-redshift SNe la sample. Figure 3 
distills the t^xp differences taken from Figure 2 into equivalent Gaussian standard deviations for 
the difference in texp between the high-redshift and low-redshift SNe la samples. These plots 
demonstrate that for tjoin < —10, the high-redshift and low-redshift SNe la samples agree at the 
l-cr level or better. For tjoin less than ^ — 15 days, the high-redshift SNe la sample is unable to 
place meaningful constraints on texp- 

The rise-time value quoted in Riess et al. (1999a) of texp = —19.98 ±0.15 — compared to 
texp = —20.08 ± 0.19 (statistical) obtained from our analysis — was determined for tjoin — 

-10 

days, and is plotted in Figure 2. Even at this reference epoch the disagreement between the high- 
redshift and low-redshift SNe la samples is only 1.5-cr, not the 5.8-0" difference found by Riess et al. 
(1999a). The value of texp = —17.6 ± 0.4 days given in preliminary analysis of the high-redshift 



^SN 1997aj was excluded due the presence of several highly deviant points in its light curve (including large 
deviations within a given night) which for some combinations of texp and tjoin produced fits with greatly improved 
values of x^, but unacceptably large values of stretch. Inclusion of SN 1997aj gave longer rise times, in better 
agreement with Riess et al. (1999a), and reduced the rise-time difference by 0.9 days compared to our results in 
§2.4. Thus, although SN 1997aj was found to reinforce the findings discussed below, the most conservative choice 
was to eliminate this SN. 
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Fig. 2.— The 1-, 2-, and 3-ct (68.3%, 95.4%, and 99.73%, in progressively lighter shades of 
blue) confidence regions for the conditional probability of iex-p given tjom, -P(iexp|ijom)i for the 
high-redshift SNe la sample. The small solid circles are our results for the maximum-likelihood 
texp for a given tjom for the low-redshift sample. The uncertainties in tg^p for these points vary 
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-16. For most values of tjomi the values of 



from 0.14 days at t 

texp for the low-redshift and high-redshift samples are in fair to good agreement. The large solid 
diamond represents the best-fit texp for the low-redshift supernovae found by Riess et al. (1999a), 
for tjoin = — 10 days. The large solid square is the best-fit texp for the high-redshift supernovae 
from the preliminary work of Groom (1998) and Goldhaber (1998). The best-fit values found in 
this work for tjoin = — 10 days are in good agreement with these previous studies, especially given 
the larger uncertainties (mainly for the high-redshift SNe la sample) we find with our analysis. 
Note that the region t^xp < tjoin (shaded in black) is physically excluded by the requirement that 
a supernova light curve be a single- valued function of time. 
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Epoch # 
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-16 -14 -12 -10 -8 -6 -4 

Fig. 3. — The equivalent number of Gaussian standard deviations by which the best-fit texp for the 
low-redshift SNe la sample differs from that for the high-redshift SNe la sample, for each value of 
tjoin- This calculation accounts for the uncertainties in texp from both datasets. The difference of 
1.5a at tjoin = — 10 days is shown as a point of reference, as the studies of Ricss ct al. (1999a); 
Groom (1998); Goldhabcr (1998) have fit laws up through -10 days (to data that was pre-aligned 
by light-curve epoch and stretch). Note that for any value of tjoin the differences are much less 
than the 5.8 standard deviations found by Riess et al. (1999a). 
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sample by Groom (1998) is also plotted in Figure 2. As Figure 2 shows, the main difference between 
our finding and that of Riess et al. (1999a) lies in different best-fit values and larger uncertainties 
for the high-redshift SNe la sample (differing by —0.7 days at tjoin = —10 days). The uncertainties 
are larger, especially for the high-redshift SNe la sample, when uncertainties in the light-curve fit 
parameters, fmax^ tmax^ s (and to a lesser extent amongst the photometry points) are fully taken 
into account. These larger uncertainties come about because the individual SNe la are given the 
proper freedom to adjust to templates away from the global best-fit template. Previous analyses 
have artificially suppressed this freedom, and have therefore underestimated the uncertainty on 

^exp- 

Given the large uncertainty in t^xpi potential perturbations from the systematic effects dis- 
cussed in the next section, and the fair to good agreement in tg^p between the low- and high-redshift 
SNe la for reasonable values of tjoin-, we consider a detailed analysis of the best tjoin unwarranted. 
Riess et al. (1999b) found that per degree of freedom deteriorated for their fits for tjmn > —8 
days, indicating that the simple Q;,tei;p model is not appropriate later than —8 days. A cursory 
examination of the joint probability, P{texp, tjoin), for our fits showed that the low-redshift SNe la 
sample prefers tjoin ^ —8 days, where our analysis finds a modest disagreement between the low- 
redshift and high-redshift supornovac. However, the early low-redshift SNe la observations prefer 
a slightly different tjoin', P{texp, tjoin) based on observations having t — tmax < —6 days gives a 
preferred tjoin ^ — 12 days, where high- and low-redshift rise times agree quite well. This mild ten- 
sion within the low-redshift SNe la sample with regard to the preferred tjoin is somewhat less than 
the 2-(7 level. A similar, but weaker, situation is found for the high-redshift SNe la sample. This 
is not a complete surprise; as the following section demonstrates, there are systematic variations 
in the late-time light-curve behavior of SNe la (such as SN 1994D from the low-redshift SNe la 
sample) which can affect the preferred rise time. Furthermore, a best fit value of tjoin depends 
not only on the rise-time behavior, but also the accuracy of the modified Leibundgut template for 
tjoin <t — tmax < "4 days (the latest tjoin tested). The relative probabilites at different values of 
tjoin include a contribution from the a,t1^p model for t < tjoin and from the Leibundgut template 
for t > tjoin- Because the parameters {t„iax and s) for the modified Leibundgut template are driven 
largely by points with t > —4 days, any early-time mismatch between the modified Leibundgut 
template and the data will degrade the quality of the fit a different amount for different values of 
tjoin- This effect should only be of importance for tjoin later than about —10 days, where the data 
are better and where the best-fit a, t^^p curves begin to depart from the (full) modified Leibundgut 
template. As things stand, the goodness of fit changes imperceptibly with tjoin for the high-redshift 
SNe la sample. 



3. Systematic Effects 



Given these findings from the statistical analysis it is clear that there is a reasonable consistency 
between the rise times of the high- and low-redshift SNe la. However, it is important to explore the 
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possibility of systematic effects which have the potential to drive a fit to another location and/or 
increase the error bars further. One such effect arises from application of the stretch relationship 
when fitting an observed light curve with a given template. 

As mentioned in §2.1, the stretch method works particularly well up to t +25 days past 
maximum. After this point the light curve of a SN la leaves the photospheric phase and enters 
into the nebular phase. This is marked by a bend in the light curve between +25 and +35 days 
after maximum light where the rapid drop from peak brightness slows down into an exponential 
decline of the light curve. Since this exponential decline is governed mostly by the radioactive 
decay of ^^Co to ^^Fe one would not expect it to "stretch" like the earlier portion of the light 
curve. In fact, as seen in Leibundgut (1988), the slopes of the declines are very similar for a wide 
range of SNe la light-curve widths. This highlights one of the current limitations of the stretch 
method; the entire template, regardless of epoch, is stretched to fit the data. This is not just a 
problem for the stretch method, but for any of the current SN la template fitting methods, which 
all employ a one-to-one correlation between peak brightness and the shape of the light curve. This 
is a small effect compared to the peak flux and the typical photometric uncertainties in current 
low- and high-redshift data sets. However, it is important to consider its effect specifically on the 
measurement of the rise time. 

The amplitude with respect to peak of the aforementioned exponential decline varies among 
SNe la. It turns out that the stretch method can compensate somewhat for these differing ampli- 
tudes, providing better fits in the sense, but at the expense of introducing a possible bias in 
s. Since the amplitude variations during the exponential decline become apparent at brightnesses 
similar to those on the rising portion of the light curve being studied here, and since the data are 
generally much better for the later portion of the light curve, the late-time light-curve behavior 
may bias determination of the rise time. The eff'ect of this bias on the template fitting method was 
studied via a Monte Carlo simulation, as described below. 

Figure 4a shows the modified Leibundgut template along with two other templates derived from 
the SNe la 1986G and 1994D (Phillips et al. 1987; Meikle et al. 1996; Patat ct al. 1996). These 
supernovae were chosen because, among those SNe la with good late-time data, they produced 
the largest deviations from the modified Leibundgut template in the tail of the light curve. To 
produce these templates for the Monte Carlo simulations the data through day +15 for SNe 1986G 
and 1994D were adjusted to fit the modified, unity-stretch, Leibundgut template. The resulting 
adjustments were then applied to the data beyond day +15 using the stretch method. These 
adjusted late-time data were fit with a smooth curve through the bend in the light curve, followed 
by an exponential decline. These late-time curves were then mated to the modified, unity-stretch, 
Leibundgut template for t < +15 days to form complete templates. 

Figure 4b shows the normalized ensemble photometric error for both the high-redshift and low- 
redshift SN la samples in 7 day bins from —35 < t — tmax < +75. This indicates how accurately 
a light curve would have been measured had all the observations come from just one supernova. 
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Fig. 4. — a) The different light-curves shapes used to study the systematic uncertainty due to 
temporal sampling and intrinsic light-curve deviations. The black curve corresponds to the stan- 
dard modified Leibundgut template. The blue curve shows the late-time deviation exhibited by 
SN 1986G, while the green curve shows the late-time deviation for SN 1994D. b) The effective 
normalized ensemble uncertainty as a function of time for the high-redshift SNe la (black) and 
low-redshift SNe la (red) samples. If the observations from each of these samples had come from 
a single high-redshift or low-redshift supernova, this would reflect how well that light curve is 
determined. Although the uncertainty is not strongly dependent on light-curve phase for the high- 
redshift SNe la, the scale over which the data can be effectively averaged — set by the steepness 
of the light curve — leads to a poorer determination of the early, rising, portion of the light curve. 
Note also that the uncertainty is modulated slightly by lunar phase, since the SNe la are typically 
discovered shortly before maximum light and one week prior to new moon, and cannot be observed 
near full moon. 
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Similarly, provided the stretch method works sufficiently well, and fmax-, tmax, and s are known, 
this would be the accuracy of a stretch-corrected composite light curve. Note that the high-redshift 
data are of consistent quality through 50 days after maximum light, which enables the high- 
redshift SNe la sample data to constrain the fit to a template over a large range in time with nearly 
equal weight. However, this makes the high-redshift SNe la data susceptible to a systematic bias on 
the rise time due to possible deviations from the stretch fitting method for t > 40 days for deviant 
light curves like those shown in Figure 4a. 

The Monte Carlo simulation performed to test for such a bias created simulated light-curve 
photometry data for the three sets of supernovae based on the templates seen in Figure 4a. Each set 
was comprised of ~ 100 different realizations of each of the supernovae in the high-redshift SNe la 
sample based on their individual temporal sampling and associated photometry errors. All of the 
generated supernovae were created with the following input parameters: s = 1.0, tjoin = —10.0 
days, texp = —20.0 days, tmax = 0-0 and fmax = l-O- The resultant light curves produced in each 
set were fit with the modified Leibundgut template, surfaces of texp and tjoin were created for 
each of the fits, and within a set these surfaces were added together to find the global minimum. 

The results of these simulations are given in Table 1. It is apparent that given a set of SN la 
observations like those available from the high-redshift SNe la sample, a fit for texp can be biased 
by 2 — 3 days in either direction if all the observed SNe la have deviant late-type light curves like 
SN 1986G or SN 1994D. To allow direct comparison with Figure 2, these same simulations were 
used to determine the best values of texp for input templates with tjoin fixed at —10.0 days. For this 
case we found texp = —19.8, —17.5, and —22.5, when the data were simulated using the Leibundgut, 
SN 1986G, and SN 1994D templates, respectively. This shows that systematic errors in texp are 
large even when tjoin 

is held fixed. While the SNe la template light curves used to simulate the 
high-redshift data can be thought of as extreme cases, at present the exact nature and frequency 
of such deviations at this light-curve phase is poorly quantified due to a lack of high-quality, well- 
sampled observations over peak and through day ~ -|-60 for nearby supernovae. Therefore, this 
result should be taken as a rough upper limit on the systematic error on texp due to temporal 
sampling and our current limited understanding of how the stretch relationship should be applied 
at late times. 



Table 1. 
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Leibundgut 
SN 1986G 
SN 1994D 


-20.00 
-20.00 
-20.00 


-10.0 
-10.0 
-10.0 


-19.76(62) 
-18.08(59) 
-23.62(19) 


-8.01(42) 
-6.81(44) 
-8.84(20) 


-F0.006(011) 
-F0.039(012) 
-0.020(009) 
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4. Cosmological Implications 

Assuming all SNe la have rise times similar to that found by Riess et al. (1999b) from good 
early-time photometry, a light-curve template with texp ~ — 20 days and tjoin ~ — 10 days might be 
a better template for use in fitting the light curves of SNe la at all redshifts. This raises the question 
of whether such a change from the modified Leibundgiit template to a Riess-like template would 
alter the corrected peak magnitudes determined in Pcrlmutter et al. (1999). In comparing our fits 
to the high-redshift SNe la sample using these two alternative templates, wc find no measurable 
change in the ensemble mean corrected peak magnitudes. We also find that no individual SN la 
changed by more than 0.02 magnitudes. 

Another obvious question, addressed by the simulations of §3, is whether systematic variations 
in late-time light-curve behavior can affect the cosmological results of Perlmutter et al. (1999). 
In the last column of Table 1 we list AMg'^'', the change in the ensemble stretch-corrected peak 
magnitude for each dataset determined using the stretch-luminosity relation of Perlmutter et al. 
(1999). These changes (—0.039 < AM^"' < 0.020) arc small, and less than the systematic biases 
already considered in Perlmutter et al. (1999) (0.05 mag). Given the fact that these simulations 
represent the most extreme deviations encountered with our fitting method, we conclude that this 
bias has no effect on the determination of the cosmological parameters from SNe la. 

5. Conclusions &; Discussion 

We find no compelling statistical evidence for a rise-time difference between nearby and distant 
SNe la, and therefore no evidence for evolution of SN la. We do find that for the high-redshift 
SNe la sample, temporal sampling coupled with real deviations of SNe la light curves at late- 
times could systematically bias the inferred rise time by 2 — 3 days. Even if present, these biases 
cannot dim the peak magnitudes by more than 0.02 magnitudes nor brighten them by more than 
0.04 magnitudes even in the extreme cases that all the distant SNe la have late-time light curves 
like SN 1994D or SN 1986G, respectively. This leaves the cosmological results of Perlmutter et al. 
(1999) unchanged. Due to the large statistical uncertainties and possible systematic effects, we 
conclude that the extant photometry of high-redshift SNe la are in fact poorly suited for placing 
meaningful constraints on SN la evolution from their rise times. 

If future studies using better early-epoch data (such as that expected from the SNAP satel- 
lite^) were to find significant rise-time differences between nearby and distant SNe la, would this 
invalidate the use of SNe la as calibrated standard candles? This is a very complicated question. 
However, at least some models suggest that variations in the early rise-time behavior may be very 
sensitive to the spatial distribution of ^^Ni immediately after the explosion. Such differences would 



^See http://snap.lbl.gov for information pertaining to the SuperNova Acceleration Probe. 
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diminish as the SN la expands and the photosphere recedes, meaning that rise-time variations 
wouldn't necessarily translate into differences in peak brightness (Pinto 1999, private communica- 
tion). Careful measurement of the rise time and the peak spectral energy distribution of individual 
SNe la will have to be carried out to address this question (see Nugent et al. (1995a,b) for a full 
description of the interplay between the rise time and the spectral energy distribution on the peak 
brightness of a SN la). It may even prove possible to use the rise time as an additional parameter 
to improve the standardization of SNe la. 

We close with some general observations concerning the issue of SN la evolution. The peak 
brightnesses of SNe la are determined at some level by the underlying physical parameters of metal- 
licity and progenitor mass, whose mean values can be expected to evolve with redshift. Nonetheless, 
there should exist nearby analogs for most distant SNe la since there is active star formation and 
a wide range of metallicities within nearby galaxies (Henry and Worthey 1999; Kobulnicky and 
Zaritsky 1999). The existing empirical relations between intrinsic luminosity and light-curve shape 
are able to homogenize almost all nearby SNe la. This implies that SNe la with some finite (but 
as yet poorly quantified) range of metallicities and progenitor masses can be used as calibrated 
standard candles. This forms the basis for using SNc la at high-rcdshift to probe the cosmology. If 
there is a dominant population of SNe la whose members are undcrluminous for their light curve 
shape at z ~ 0.5, as would be required to explain current observations in terms of evolution, there 
should be nearby examples of these SNe la. Such SNe la are not predominant among nearby SNe la, 
as almost all nearby SNe la obey a width-brightness relation. For such SNe la to predominate at 
z ~ 0.5 while being rare nearby requires a large reduction in their rate. Searches for SNc la con- 
ducted using exactly the same CCD-based wide-area blind-search methods used by the SCP find 
that the SNc la rate per comoving volume element does not change significantly between z < 0.1 
(Aldering 2000), z ~ 0.5 (Pain et al. 1996, 2000, in preparation), and z ~ 1.2 (Aldering et al. 2000, 
in preparation). For the global rates to stay roughly constant while the rate of such hypothetical 
subluminous SNe la changes by an order of magnitude would be remarkable. For instance, a shift 
from Pop II progenitors at z ~ 0.5 to Pop I progenitors nearby would result in suppressed rates 
at z ~ 0.5. This is due to the fact that Pop II stars are a minor contributor to the luminosity 
density out to z ~ 0.5 (Shimasaku and Fukugita 1998). Quantifying these arguments is beyond the 
scope of this paper, so we do not claim they as yet place a bound on SN la evolution. However, 
such arguments should be borne in mind when weighing the likelihood that the calibrated peak 
brightnesses of SNe la evolve. These arguments can also provide a partial basis for rigorous testing 
of the SN la evolution hypothesis. 
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